pub fn new_birthday_probability(n: u32) -> f64 {
    let year_days = 365.;
    let mut possibility = 1.;
    for k in 0..n {
        let p = k as f64;
        possibility *= (year_days - p) / year_days;
    }
    1. - possibility
}
